Control of quantum dynamics by optimized measurements 
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Abstract 

Quantum measurements are considered for optimal control of quantum dynamics with instan- 
taneous and continuous observations utilized to manipulate population transfer. With an optimal 
set of measurements, the highest yield in a two-level system can be obtained. The analytical solu- 
tion is given for the problem of population transfer by measurement-assisted coherent control in a 
three-level system with a dynamical symmetry. The anti-Zeno effect is recovered in the controlled 
processes. The demonstrations in the paper show that suitable observations can be powerful tools 
in the manipulation of quantum dynamics. 
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I. INTRODUCTION 

Control of quantum processes has attracted considerable attention, both theoretically Q, 
and experimentally |7|, |8|, |9J . Most studies of quantum control are concerned 
with shaping a laser pulse to achieve an expected goal. However, a coherent laser pulse is not 
the only driving force for quantum dynamics. Incoherent driving forces, such as laser noise, 
decoherence from the environment, and quantum observations, can also influence quantum 
dynamics. A natural general expectation is that the later influence will be deleterious 



toward achieving control [10j. However, recent studies [111. 112] have shown that controlled 
quantum dynamics can survive intense field noise and decoherence, as well as even cooperate 
with them under special circumstances [yj]. Under these special conditions, it is possible 
to meet the target goal more effectively with the help of laser noise and environmental 
decoherence. Incoherent non-unitary control of quantum systems by a suitably optimized 
environment (e.g., incoherent radiation, a gas or solvent, a cloud of electrons, atoms or 
molecules, etc.) was suggested as a supplement to coherent control to provide a general tool 
for selective manipulation of both the Hamiltonian and dissipative aspects of the system 

dynamics HQ. 

Both the outcome and back-action from quantum measurements could be used to control 



quantum processes. In standard closed- loop optimal control [16|, the quantum system is 
non-selectively measured when the desired evolution ends, and the outcomes of the mea- 
surements are used by a learning algorithm to optimize the laser pulse. Measurements were 
also used to map an unknown mixed state onto a known target pure state 17|] . Some inves- 
tigations considered exploiting the back action from the quantum observations, and control 



by measurement plus dynamical evolution was proposed 



181 ] . The control of the population 



branching ratio between two degenerate states by continuous measurements was treated 19], 

21]. 
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and the effect of nonoptimized measurements on control by lasers was investigated 

Numerical simulations have been performed to investi gate observations serving as indi- 
rect controls in the manipulation of quantum dynamics [22j. Optimal control fields were 
shown to be capable of cooperating or fighting with observations to achieve a good yield, 
and the nature of the observations may be optimized to more effectively control the quan- 
tum dynamics. Quantum observations also can break dynamical symmetries to increase the 
controllability of a quantum system. The quantum Zeno and anti-Zeno effects induced by 
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observations are the key operating principles in these processes. When instantaneous obser- 
vations are the only forces to drive population transfer in a two-level system, the maximal 
population transfer induced by a finite number of measurements has been found, and the 
quantum anti-Zeno effect is recovered in the limit of infinitely frequent measurements 2jJ. 

In this paper, we further explore the utility of quantum measurements as controls in the 
manipulation of quantum dynamics. Here we assume that any projection operator may be 
observed in a suitably performed experiment. Analytical solutions and upper bounds of 
several controlled processes are found. The remainder of the paper is organized as follows: 
Section II reviews the main concepts of instantaneous and continuous measurements, which 
are used as controls in this paper. The analytical solutions for instantaneous and continuous 
observations in a two-level system are explored in Sec. Ill and IV, respectively. The maximal 
measurement-assisted population transfer in a system with dynamical symmetry is obtained 
in section V. A brief summary is presented in section VI. 

II. QUANTUM OBSERVATIONS 

Quantum measurement serves as an incoherent driving force, and there are two general 
types of quantum measurements: instantaneous von Neumann measurements and continuous 
measurements. A characteristic feature of quantum systems is that their measurement 
unavoidably affects the associated dynamics. The well known manifestation of this back 



reaction is the uncertainty principle |24|. The influence of a measurement is revealed in a 
more direct way through a change of the measured system state. In von Neumann axiomatic 
quantum mechanics it is postulated that any measurement gives rise to an abrupt change 
of the state of the system (instantaneous measurements) under consideration and projects 



it onto an eigenstate of the measured observable [25|. The measurement process follows 
irreversible dynamics, and causes the disappearance of coherence of the system's state: 
the off-diagonal elements of the density matrix decay, or the phases of the wavefunction 
amplitudes are randomized. Density matrices are adopted to describe the states of controlled 
systems, because the nonselective quantum measurements in the paper are applied to an 
ensemble of quantum systems. 
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A. Instantaneous measurements 



An ideal measurement occurs at one instant of time or a sequence of such observations 



may be performed at different instants 



251 ] . Let Q = YliiQiPi ^ e an observable with qi 



being the eigenvalue of Pi, which is a projector such that P%Pj = Pi^ij and £\ Pi = 1. The 
measurement of Q converts the state p of the system just before the measurement into the 
state 

k 

A projection operator P satisfies P = P 2 , and its spectral decomposition may be written 
as P = qoPo + q\P\ with the two eigenvalues being q = and q% — 1 and two corresponding 
projectors being Pq = 1 — P and P\ = P. Thus, according to Eq. (TjQ), observation of the 
operator Q = P transforms the density matrix p of the system into p' given by 

p' = p P (p) = P pP + P lP P 1 (2a) 
= p-[P, [P, p}} , (2b) 

so [P, [P, p]} is the "kick" by an instantaneous observation of the projection operator P. 
B. Continuous measurements 

There are two equivalent theoretical formulations of continuous quantum measure- 
ments 26|. One of them is based on restricted path integrals (RPI) and the other one 



on master equations (ME). For simplicity, we adopt the latter formulation. Corresponding 
to a continuous measurement of a single observable, the master equation has the form j^ : 

p=-i[H,p]-± K [A,[A,p]\. (3) 

Here, H is the Hamiltonian of the measured system, A is the measured operator. Equa- 
tion ([3]) is similar with the equation describing a system interacting with the environment, 
so we could call k as the "strength" of the observation. 



C. Quantum ZENO and Anti-ZENO Effect 

Prevention of a quantum system's time evolution by means of repetitive, frequent observa- 
tions or continuous observations of the system's state is the quantum Zeno effect (QZE). The 



QZE was proposed by Misra and Sudarshan 28] and was experimentally demonstrated 29] 



in a repeatedly measured two-level system undergoing Rabi oscillations. A time- dependent 
observable projection operator inducing up to 100% transfer from one state to another 
state {30] is the quantum anti-Zeno effect (QAZE). The impacts of QZE and QAZE opera- 
tions are the key processes explored in this paper to help control quantum dynamics. 

III. TWO-LEVEL SYSTEM CONTROLLED BY INSTANTANEOUS MEASURE- 
MENTS 

This section presents the analytical solution for the control of population transfer in a 
two-level system by optimized instantaneous measurements. The objective is to drive the 
population from the initial state po = |0) (0| to the final state pf = |1) (1|. An observable 
Q has the form Q = qiP^ + q%P\ [ , where q± and q2 its the eigenvalues, and P^ and 
P^ are the corresponding projectors. From Eq. (T5]), it's easy to see that the measurement 
of Q is equal to the measurement of the projector P[ , since P% — 1 — P\ ■ Thus, it is 
sufficient to consider the measurement of projection operators in this case. 

A sequence of N instantaneous projection observations, specified by the operators 

Pk = M , (4a) 
|^) = cos^|0) + e^sin^|l). (4b) 



are performed at times J 1 *;, k — 1, • • • , N. Parameters and 9k in Eq. (14b j) are limited to 
the range 

"2 £ f < 2" (5a) 

< 9 k < 7T, (5b) 

since the projection operator Pk does not depend on the phase of \ipk)- The operators Pk 
parameterized by a& and 9k, k = 1, • • • , N, are optimized to maximize the yield, 

Y N [P lt --- ,P N ] = (l\p N \l). (6) 

Here pj^ is the density matrix after performance of N observations given by the iterative 
equation 

Pk = Pk-i ~ [ p k, [Pk,Pk-i]) , k — ■ ,N (7) 



with p = |0) (0| and Pk described in Eq. ( I4al) . We have neglected the free evolution between 
measurements, which is easy to include via a transformation between the Schrodinger picture 
and interaction picture, 

p k = e m « T *pi S) e- m « T \ (8a) 
Pk = e ^T kp (s) e - lH0 T K (gb) 

Here p k is the density matrix in the Schrodinger picture governed by the iterative equation, 

p (S) = e --o(T fc -T fc _ l){p (5) _ [p (*) ? [p (5) j p f_) l]]}e -o(T fc -T fc . l); k = lt ... iNi (9) 

which includes the free evolution. 

The density matrix p^ is Hermitian with unit trace. Hence, it can be expressed in the 

for 

^f 1 -***). do) 

\ Z N Y N j 

It is easy to establish the following solution to the iteration in Eq. (JTj) 

Y k = (l\p k \l) = i(l-cosaiCi 2 C 2 3---C fc _i ifc cosa /c ), (11a) 
Z k = (l|pfe|0) = ^e l6k cosaiCi 2 C23 • • • C fc _i ifc sin (lib) 

with the coefficients C mn given by 

C mn = cos a m cos a n + cos (9 m - 6 n ) sin a m sin a n . (12) 

Therefore, 

Ijv = - (1 - cosaiCi 2 C 2 3 • • - Cat- i, v cos a N ) (13) 

is the yield from iV observations. 

We will now determine the maximum value of Yn, which is a function of variables 9 k and 
aifc. The inequality ( l5bl) yields — 7T < 9 m — 9 n < tv and 1 > cos (0 m — 9 n ) > —1. Setting 
to zero the derivative of Yjv with respect to ^ gives sin {Ok ~ #fc-i) = for = 2, • • • , N . 
Hence Y^ reaches its maximum when 

9 1 = 9 2 = --- = 9 N (14) 
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Therefore, after optimization with respect to 6k, Yn can be written as an function of a^, 



k = 1, 



,N, 



= — [1 — COS OL\ COS («i — OJ2) • • • COS (cfcjv-l — Q!jv) COS«Ar] 



niV 
fe=0 



COSy^ 



(15a) 
(15b) 



with ip = n — «!, <fi — q;i — a 2 , • • ■ , y^-i = a iv-i — a iv and yjjv = ajv- It is easy to verify 
that the second derivative of the function / (2) = In cos (x) is negative, so it is a concave 
function. The inequality, 

nL 008 ^- ( cos ivTi) ' (l6) 

can be established by the majorization inequality [3l|] for concave functions, 



V M f(Xl.) Xh 



M 



M 



(17) 



Hence, reaches its maximum value 



Y, 



(O) 



iV 



7T 



COS • 



A + 1 



iV+l 



(18) 



when <£> = ■ ■ • = <£n = ]y+T- The solutions are consistent with what was found in [23(. The 
QAZE is recovered in the limit of an infinite number of observations, 



lim Y, 



(o) 



N 



1. 



(19) 



IV. TWO-LEVEL SYSTEM CONTROLLED BY CONTINUOUS MEASURE- 
MENTS 

In this section, the quantum dynamics of a two-level system is controlled by suitable 
continuous measurements. Here we assume that it is possible to continuously measure any 
time-dependent projection operator P(t). In the interaction picture, the dynamics of the 
continuous observation process is described by 

p (t) = ->y£ (t) P (t) = -7 [P (t) , [P (t) , P («)]] , (20) 
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where C (t) is a super-operator acting on the density matrix p(t), and 7 is the constant 
strength of the observation. The projection operator P (t) is specified by 

P(t) = \^(t))(^(t)\, (21a) 
IV, (*)) = cos ^ |0) + e ie « sin ^ |1) , (21b) 

where a (£) , 6* (t) are functions of time t to be determined. The goal is to optimize the 
objective functional, as the yield at final time Tf, 

Y (T f ) [P (t)] = Y (Tf) [a (t) , 6 (t)] = (l\p (T f ) |1> , (22) 

where the system is initially populated on state |0). 

Eq. (l20|) appears insoluble for general functions of a (t) and 8 (t). First consider only a 
simple case, with zero phase and a (t) taken as linear in time, 

a (t) = A— + B (23a) 
6 (t) = 0. (23b) 

The final yield of this case may be explicitly worked out as 

Y (T f ) = -- -e^' { cos A cosh 5 + [7' cos (2B + A) + A sin A] j , 

2 2 1 J 

where 7' and 5 are dimensionless parameters defined by 



(24) 



7' = l/T f (25a) 

5 = vV 2 - A 2 . (25b) 

Eq. (1241) reaches its maximum value when 

2B m + A m = n (26a) 

7' sin A m = A m , (26b) 

in which case the optimal value of the yield is 

(Tf) = i (1 - e -7'(i+cosA m ) cosAm j . (27) 
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Figure 1 depicts the variation of the optimal A m and B m with respect to 7', and we conclude 
that 



A m = 0, for 7' < 1 (28a) 
A m -> 7r, for 7' -> 00. (28b) 

It follows from Eq. fT2~T|) that a complete population transfer is attained when 7' in 
Eq. (125aj) approaches infinity. Thus, from Eq. (125a[) . increasing the observation strength 
7 and the final time are equally effective in enhancing the control process, and the QAZE 
is recovered in the limit of infinite observation strength, or final time. 

We now assess whether the linear solution in Eq. ( |23i) is optimal with respect to all 
possible forms of a(t) and 0(t). To verify that this is the case, we start from Eq. (|2"U|) 
and consider the variation of p(t) with respect to a (t) and 9 (t). The general variation of 
Eq. (|2DP gives 

^£^1 = -7 [C it) 5p (t) + 5C (t) p (t)} . (29) 

and for a driving variation 5a (t), it is easy to verify that the solution of the above equation 
is 

5p (t) = —7 / [p a (t, t)] 5a (r) dr, (30a) 
Jo 

p a (t,r)=U(t,r)[^-p(r)}, (30b) 



—7 / C (y) dv 



where U (T, t) is a time-ordered exponential 

U (t, t) = exp + 

Hence p a (t, r) is the solution of the differential equation 

dp a (t, t) 
dt 

with the initial condition 



(31) 



-jC(t)p a (t,T) (32) 



dC (t) 

Pa (t, t) = ^ p (r) . (33) 

It is evident that p a (t, r) is real, symmetric and traceless under the assumption of Eq. ff23l . 
hence we can set 

' -Y a (t,r) Z a (t,r) 

Pa (t,T)=\ • (34) 



It is easy to yield the result 



Y a (t, r) = ~ sin A ^1 - ^- ) exp 



A(t- 2r) cot A - t esc A 



Tf 

So the variation of the final yield with respect to a (t) is zero, 



(35) 



gy(.o) ( Tf j = _ 7 f 1 Ya (r f , t) 5a (r) dt = 0. (36) 
Jo 

Using the same procedures taken above, we can prove that the variation of the final yield 
with respect to the phase function 9 (t) is also zero. Hence, the linear solution is an optimal 
solution. 

In addition to the analysis above, we performed numerical simulations, where the goal was 
optimization of the yield Y by means of an evolutionary algorithm approach to determine 
a(t) and 0(t). The optimization procedure was conducted freely, without any preliminary 
assumptions on a(t) or 9(t), and without any constraints on their values during the search. 
For this purpose, we applied the covariance matrix adaptation evolution strategy (CMA- 



ES) |32|, |33j to the task. 



optimization problems 



ie latter algorithm is very efficient for treating continuous global 



351 ] . It has been successful for handling correlations among object 



variables. Fig. 2 depicts the best yield from performing full optimization with respect to 
a (t) and 9 (t), when the assumptions in Eq. ( 1231) are not applied. It can be concluded that 
the solution is globally optimal. 

The best yield from continuous measurements, Eq. (1271) . seems very different to the best 
yield from instantaneous measurements, Eq. ffl8|) . However, their asymptotic forms, 

2 

4 0) ~l-^> N^^, (37a) 
Y<°^T f ) ~ 1 - (37b) 

are very similar. Hence the best yield from continuous measurements with measurement 
strength 7 is very close to the best yield from a sequence of iV ~ jTf/2 instantaneous 
measurements, when 7TJ ^> 1. 
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V. OPTIMAL POPULATION TRANSFER IN A SYSTEM WITH DYNAMICAL 
SYMMETRY BY MEASUREMENT-ASSISTED COHERENT CONTROL 



In this section we consider a system whose free Hamiltonian H and dipole moment ji 
are given by 



( 



1 
2 
3 



\ 



/ 



V 



1 

1 1 
1 



\ 



/ 



The system is initially prepared in the ground state j^o) = |0) at t — 0. The control goal 
is to transfer as much as possible of the population from the ground state |0) to the first 
excited state |1) at the target time T > using as controls a coherent electromagnetic field 
e(t) during the time period [0, T] and a single measurement of the projector P — |0)(0| (or 
P 2 = |2)(2|) at a time h e (0,T). 



The symmetry of the system implies (see Ref. 36] ) that the coefficients of a pure system 
state \ip t ) = Co(t)\0) + C\(t)\l) + C2(t)\2) satisfy the following relation upon evolution under 
only the action of a coherent field, 



C (t)C 2 (t) 



Cf(t) 



C (0)C 2 (0) 



C?(0) 



for any t > 0. 



(38) 



If the initial state is \ipo) = |0), then Ci(0) = 6*2(0) = and Eq. fl38|) becomes the following 
relation for the coefficients Ci(t): 



Cl{t) = 2C (t)C 2 (t) , for any t > 0. 



(39) 



This relation was used in 22j to conclude that transferring more than 50% of the population 
from the level |0) to the level |1) is impossible using only a coherent control field. 

Measurements performed on the system can break the dynamical symmetry thereby 
allowing for exceeding the above 50% population transfer limitation. Numerically, the 
measurement-assisted control problem for this system was investigated in Ref. 22j, where 
transferring 66.9% of the population to the level |1) was obtained with a coherent control 
field assisted by a single measurement of Pq. In this section we analytically treat this control 
problem to find the upper bound on the maximal population transfer to the level |1), which 
is found to be approximately 68.7%. We also explicitly find the Rabi Frequencies of the 
optimal pulses, thus providing a complete analytical solution to this problem. 
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The control process consists of the following three steps. First, the system evolves under 
the action of a coherent field during the time interval [0, ti). Second, at the time t — t\ a 
non-selective measurement of Pq is performed on the system, which transforms the system 
state in accordance with the von Neumann scheme. Third, the system evolves during the 
time interval [ti,T] again only under the action of a coherent field. 

Spontaneous emission during the first and third steps is neglected in this consideration. 
Therefore the system's dynamics under the action of an electromagnetic coherent field during 
the first and third steps can be described by optical Bloch's equation without relaxation 
terms: 

^ = -i[H,p(t)}. (40) 

Here the Hamiltonian H = Q(t)\0)(l\ +Q(t)\l)(2\ +h.c. is determined by the Rabi frequency 
Q(t) of the electromagnetic field e(t). The symmetry of the system implies that the Rabi 
frequencies for the transitions |0) «-> |1) and |1) «-> |2) are the same. 

Without loss of generality, it is sufficient to consider, constant Rabi frequencies during 
each step of the control. If the Rabi frequency Q is time independent, then the Hamiltonian 
has the form H = H(Q) = fi|0)(l| + fi|l)(2| + h.c. and the solution of (gQJ) with the 
initial condition p(t ) = p is pit) = U(t - t )p W(t - t ), where U{r) = e ~ irH{n) . The 
Hamiltonian H(Cl) can be written in terms of the vector \Q) = [O|0) + fi*|2)]/\/2|f2| as 
H(Q) = y/2\n\[\l)(n\ + (here tt* is the complex conjugate of fi). One has 

[H(n)} 2 = 2|fi| 2 [|i)(i| + |n><n|] ^ f [H(n)} 2n = (V2|fi|mi>(i| + |n><n|] 
[H(n)f = 2\n\ 2 H [ [H(n)] 2n+1 = {V2\n\) 2n+1 [\i){n\ + 

This gives 

U(t) =P& + cos(V2|0|t) (jl)(l| + -ism(V2\n\rj (\l)(tt\ + (41) 

where P^ = I— |1)(1| — |^)(^| = |^)(^| is the projector onto the subspace generated by the 
vector = (fi|0) -W\2))/V2\n\. 

In the first stage of control, the initial density matrix po is transformed into p\ = U\PqU\ = 
\ipx) where j^i) = U\\0) and U\ = exp[— itiH(D,x)] is the evolution operator induced by 
the control field with some Rabi frequency Qi = l^ile"^ 1 . Direct calculations give = 
C |0) + C 1 |l) + C 2 |2) with 

cos(V2|^i|ti) +1 isin(v^|fii|*i) „ _ cos( v / 2|fii|ti) -1 

2 ' v/2 '2 2 
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Measuring the projector Po at the time t = t\ transforms the pure state p\ into the 
density matrix 



P2 = fiP (Pl) = PoPlPo + (I - ^0)P1(I " Po) 

= \C \ 2 P + \C 1 \ 2 P 1 + \C 2 \ 2 P 2 + C l C*\l)(2\ + C*C 2 |2)(1| 

If C\ 7^ 0, then the state p 2 is mixed and the measurement destroys the coherence between 
the levels |0) and |1) while preserving the coherence between |1) and |2). 

After the measurement the system density matrix evolves under action of a coherent field 
with some Rabi frequency Vt 2 = |^ 2 |e"^ 2 into p3 = U 2 p 2 U\. Here U 2 = exp[— it 2 H(Vt 2 )} is 
the evolution operator induced by the coherent field and £ 2 = T — t\. The density matrix 
P3 = U 2 [pp {UiPoUfy\ Jj\ can be computed using (J3TJ. The computation gives the following 
population P = ( 1 1 P3 1 1 ) of level |1) at the target time t = T, 

P = — < 5 — cos(si) — [1 + 3 cos(a;i)] cos(x 2 ) + 2 [2 sin(xi/2) — sin(xi)] sin(a; 2 ) cos(/02 — ipi) \, 

16 L J 

(42) 

where x\ = 2v / 2|^i|^i and x 2 = 2v / 2|f2 2 |£ 2 . This function is maximized by 



x* = ± 



(43a) 



x* 2 = =Farctan [ -=^- — ——] (43b) 

V6-1 J 

and by ipi, ip 2 such that ip 2 — V'i = 27rfc, k = 0, ±1, ±2, .... The maximal value is 

P max = max P = 4 • 1(T 3 ^393 -48^ + 138 + 7V&j w 68.7%. (44) 

This maximal population transfer to the level |1) can be obtained by applying a coherent 
field with Rabi frequency = x\j {2\f2t\)e 1 ^ 1 to the system during the time interval [0, £1], 
then measuring the projector Po at time ti, and finally applying a coherent field with Rabi 
frequency f2 2 = xl/(2^/2t 2 )e 1 ^ 2 during the time interval [ti,T], where ip 2 — ipi + 2nk, 
A: = 0, ±1, ±2, 

Another simple way to compute the maximal yield is with the well-known Euler decompo- 
sition of the SU{2) Lie group [371] . Suppose the system' propagation consists of the following 
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three steps: 

Pi = U}poUx (45a) 

p 2 = Pi-[P k ,[Pk,Pi]] (45b) 

Ps = U\p 2 U 2 (45c) 

where Euler's decomposition of the unitary propagators is as follows: 

U k = exp (ia k H ) exp ( i-^=H J exp {ib k H ) , k = 1, 2, (46) 



2\/2 

and 01,2, &i,2 and are six independent variables to be optimized with. Simple computation 
yields the population of the level |1) as 

(P3)n = iq ~~ cosx 2 — cosx! (1 + 3 cosx 2 ) + 2 cos (a 2 + hi) ^sinxi — 2 sin — J sin x 2 , 

(47) 

with which it is easy to derive the same maximal population transfer shown in Eq. (|4*4|) . 

VI. CONCLUSION 

This paper discusses the use of both instantaneous and continuous observations in the 
manipulation of quantum dynamics. The measurements can be viewed as direct controls. 
Two-level systems and a special three-level system are treated analytically. Solutions and 
upper bounds for the controlled processes are obtained, and they agree very well with previ- 
ous numerical simulations, and QAZE is recovered. The results are proper for instantaneous 
observations performed any number of times and continuous observations performed with 
any strength. The performance of optimal observations hopefully will become routine with 
advancing technology, as observations can be powerful tools in the control of quantum dy- 
namics. 
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Fig. 1 

FIG. 1: (Color online) Optimal coefficients A m and B m of the linear function a (t) in the projection 
operators (ref. Eq. (|21ap ) which are continuously measured to control the quantum dynamics of 
a two-level system. In the optimization process, the functions a (t) and 6 (t) in the projection 
operators are assumed to be linear and zero, respectively, and the optimal solutions are proved to 
be globally optimal. 7' (dimensionless) is multiplication of observation strength and observation 
time (ref. Eq. (125a]) ). 
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FIG. 2: (Color online) The yield as a function of the observation strength for the linear solution 
(solid line), with the best attained yield with an evolutionary search using the CMA-ES algorithm 
(squares) as a reference. The latter non-linear solutions did not exceed the linear solution's global 
optimal yield. 7' (dimensionless) is multiplication of observation strength and observation time 
(ref. Eq. (|25a|>). 
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